Generation of wakefields by whistlers in spin quantum magnetoplasmas 
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The excitation of electrostatic wakefields in a magnetized spin quantum plasma by the classical 
as well as the spin-induced ponderomotive force (CPF and SPF, respectively) due to whistler waves 
is reported. The nonlinear dynamics of the whistlers and the wakefields is shown to be governed by 
a coupled set of nonlinear Schrodinger (NLS) and driven Boussinesq-like equations. It is found that 
the quantum force associated with the Bohm potential introduces two characteristic length scales, 
which lead to the excitation of multiple wakefields in a strongly magnetized dense plasma (with a 
typical magnetic field strength Bo > 10 9 T and particle density no > 10 36 m -3 ), where the SPF 
strongly dominates over the CPF. In other regimes, namely B < 10 s T and n < 10 35 m~ 3 , where 
the SPF is comparable to the CPF, a plasma wakefield can also be excited self-consistently with one 
characteristic length scale. Numerical results reveal that the wakefield amplitude is enhanced by the 
quantum tunneling effect, however it is lowered by the external magnetic field. Under appropriate 
conditions, the wakefields can maintain high coherence over multiple plasma wavelengths and thereby 
accelerate electrons to extremely high energies. The results could be useful for particle acceleration 
at short scales, i.e. at nano- and micrometer scales, in magnetized dense plasmas where the driver 
is the whistler wave instead of a laser or a particle beam. 

PACS numbers: 52.35. Hr, 52.35. Mw, 52.25.Xz 



I. INTRODUCTION 

Lately, significant progress has been made in the field 
of collective plasma accelerators for attaining high elec- 
tron energies [T|-fTTj| . As is wellknown, an intense elec- 
tromagnetic (EM) pulse can create a wake of plasma os- 
cillations under the action of ponderomotive force fllj j . 
Electrons are then trapped into the wake and can thus 
be accelerated to extremely high energies, as high as giga- 
to teraelectronvolt scales. The idea that was presented 
by Tajima and Dawson [llj three decades ago, has now 
become a reality through their experimental verifications 
[Hldil]. Using intense laser beams or bunches of relativis- 
tic electrons have been used to excite plasma waves and 
to produce high electric field strengths (10-100 GV/m), 
opening up the possibility for compact GeV particle ac- 
celerators (see, e.g., [l| for a discussion). Recently, an 
alternative scheme has been proposed for accelerating 
electrons up to the TeV regime, using proton bunches 
for driving plasma- wakefield accelerators [l|. 

A question that naturally arises is whether or not 
there exist natural (high-energy) systems, such as astro- 
physical environments, in which wakefield acceleration 
can take place. It is certainly true that neither high- 
intensity lasers, ultra-short photon bunches, or particle 
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beams, used to excite wakefields in the laboratory regime, 
are readily available in the astrophysical settings, espe- 
cially in dense plasmas characterizing, e.g., the interior 
of white dwarfs and Jovian planets. In contrast to the 
above examples of typical drivers, that may exist inde- 
pendent of a plasma medium, EM whistler waves only 
exist in a plasma environment. In this manner, Chen 
et al. fT?| demonstrated the possibility of exciting large 
amplitude plasma wakefields by plasma magneto-waves, 
abundant in astrophysical settings. Later, particle-in-cell 
simulation has shown that this mechanism could well be 
valid for celestial acceleration However, no attempt 
has been made to generate wakefields self-consistently 
by whistler waves in strongly magnetized dense quantum 
plasmas with or without intrinsic spin of electrons. 

On the other hand, dense electron plasmas are degen- 
erate, and one must therefore take into account Fermi- 
Dirac statistics as well as tunneling effects [ll|. Further- 
more, apart from the statistical effects introduced by the 
spin- 1/2 nature of the electrons, the spin also generates a 
"dynamical" contribution, in the form of a pressure-like 
spin force and a spin magnetization current [171 ITcij . The 
latter can even be larger than the classical free current 
when low-frequency (If) longitudinal perturbations are 
driven by a spin-induced ponderomotive nonlinearity in 
the propagation of a short EM pulse [l9| . Naturally, the 
inclusion of the quantum statistical pressure, the quan- 
tum force associated with the Bohm potential as well as 
the spin force due to finite magnetic moment of electrons 
in the dynamical equations, will significantly change the 
collective behavior of the electrons. For example, the 
Bohm potential, arising due to the particle's wave- like na- 
ture, provides higher-order dispersion as well as the possi- 
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bility of short- wavelength (compared to the electron skin 
depth) wakefield. The spin magnetization contributes to 
linear dispersion of circularly polarized modes [13| and 
the magnetic dipole force results in a spin contribution 
to the ponderomotive force [l9j . 

In this paper, we present a theoretical investigation 
of the excitation of multiple wakefields driven by both 
the classical fCP F) [2l| as well as the spin ponderomo- 
tive force (SPF) [19| of EM whistlers in a magnetized 
quantum plasma with intrinsic magnetization. The self- 
consistent fields of the whistlers and the wake are de- 
scribed by a coupled set of nonlinear Schrodinger- and 
modified Boussinesq-like equations [22]. It is shown 
that when the SPF is comparable to CPF, the wake- 
fields can be generated in the regime of strong magnetic 
fields with Bq < 10 8 T and high-density plasmas with 
n o ^ 10 35 m~ 3 . However, multiple wakefields are also 
excited when the SPF dominates over the CPF in the 
very strongly magnetized and superdense plasmas with 
-Bo > 10 9 T, n n > 10 36 m- 3 . In the latter, the results of 
our model is that the whistler wavelength may be even 
shorter than the Compton wavelength of the electrons. 
This gives rise to a viable mechanism for wakefield ac- 
celeration in spin dominated plasmas, such as those in 
superdense astrophysical bodies, e.g., white dwarfs, neu- 
tron stars, magnetars. 



II. EVOLUTION EQUATIONS 

We consider the propagation of high-frequency (hf) 
whistler waves along an external magnetic field B = Bq2 
in a quantum plasma, taking into account the intrinsic 
spin of the electrons. The heavy ions are assumed to 
be immobile, so that the If density response occurs on 
a time-scale much shorter than the ion plasma period. 
In the modulational representation, the whistler electric 
field can be represented by E = (% — iy) E(z, t) exp(ikz — 
iujt)+c.c, where E(z,t) is the slowly varying wave enve- 
lope, ui (fc) is the whistler wave frequency (number) and 
c.c. stands for the complex conjugate. Then the basic 
equations for the evolution of whistlers consist of the fluid 
equations for the nonrelativistic evolution of spin— 1/2 
electrons given by [18|, [22( 
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and the Maxwell equations for the electromagnetic (EM) 
fields 
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Here n e , m e , v e denote the number density, mass and 
velocity of electrons respectively, B is the magnetic field 
and P e is the electron thermal pressure and S is the 
spin angular momentum with |S| = \Sq\ = fi/2, \i = 
— (g/2) /is,where g « 2.0023193 is the electron ^-factor 
and hb = eh/2m e is the Bohr magneton. The equations 
©-© are applicable even when different spin states 
(with up and down) are well represented by a macro- 
scopic average. This may, however, occur in the regimes 
of very strong magnetic fields (or a very low tempera- 
ture plasmas), where generally the electrons occupy the 
lowest energy spin states. On the other hand, for a time- 
scale larger than the spin-flip frequency, the macroscopic 
spin state (to be attenuated by a factor decreasing the 
effective value of |S| below h/2) can be well-described by 
the thermodynamic equilibrium spin configuration, and 
in this case the above fluid model can still be applied. 
However, this is not the present issue to be studied here 
rather we will focus on the regimes of strong magnetic 
fields and high density plasmas. 

The evolution equation for the whistler can be obtained 
by taking the curl of Eq. © [hence the pressure gradient 
and the quantum force in Eq. © vanish] and using Eqs. 
©-©as 
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We note that in the nonlinear interaction of hf EM waves with the If electron plasma response, the use of 
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cold plasma approximation is also justified to the fact 
that for large field intensities and moderate electron tem- 
perature, the directed speed of electrons in the hf fields is 
much larger than the random thermal speed. It can also 
be shown that the density perturbation associated with 
the hf EM wave is zero. Thus, the evolution equation 
([5]) for the whistlers do not involve contributions from 
the electron pressure and the quantum tunneling effect 
proportional to h 2 . Later, we will see that these contribu- 
tions will appear in the coupling of the If plasma response 
with the hf one through the pondcromotive force induced 
by the hf field. 

Introducing the variables B± — B x ± iB y , E± = 
E x ±iE y etc., suitable for circularly polarized waves, and 
linearizing we obtain respectively from the Faraday's law 
and the spin-evolution equation (fT3")) as [19| 
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We then linearize Eq. (J6j) and use Eq. (|Tj) to obtain 
the following linear dispersion relation for the circularly 
polarized modes [2(| 
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where ur = ck/u is the refractive index, ui p = 
g 2 | So | /4m e A 2 is the frequency due to the plasma magne- 
tization current and A e = c/uj pe is the electron skin depth 
with uj nf , 



frequency. Moreover, uj c = eB^/m e is the electron- 
cyclotron frequency and ui g — (g/2)uj c is the electron 
spin-precession frequency. We note that the frequency 
resonances occur not only at the cyclotron frequency 
(uj — > ui c ) but also due to the spin-gyration of electrons 
(uj — » ojg), although these resonances are close to each 
other. At the resonance, the transverse field associated 
with the whistlers rotates at the same speed as electrons 
gyrates around the external magnetic field. The electrons 
will then experience a continuous acceleration from the 
wave electric field. However, the detail discussion on the 
properties of the linear modes modified by the s pin mag- 
netization current can be found in the literature |20| . One 
can see that when uj c 3> uj pe , i.e., the magnetic field is 
strong enough, but ui^ is still smaller than ui pe (unless 
we consider a very high-density regime), the dispersion 
of the whistler wave are almost linear (with phase speed 
approaching the speed c of light in vacuum) over a wide 
range of wave numbers. However, it may remain non- 
linear in strongly magnetized high density plasmas de- 
pending on the parameter regimes we consider 20] . This 
nonlinear behaviors of whistlers in which E, B wave fields 
are not comparable in strength do not favor the excita- 
tion of wakefields, and we do not consider those cases 
here. 

On the other hand, in the nonlinear regime, the dy- 
namics of whistler wave envelopes can be described from 
Eq. ((6|) by the following nonlinear Schrodinger (NLS)- 
like equation [22[ • 



i (d t E + v g d z E) + (v' g /2) d 2 E — AE = 0, 
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where E = E x — iE y . v g = duj/dk is the group speed, 
d 2 ui /dk 2 is the group dispersion of whistlers and A 



\j 'riQe 2 1 'eom e denoting the electron plasma is the nonlinear frequency shift. These are given by [2^ 
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where A = 2c 2 k/u> pe + g 2 hk/4m e (uj — uj g ) , N = n e /no 
and v ez is the magnetic field aligned free electron flow 
speed. 

Note that the variables v ez and N appear due to the 



If electron plasma response, and are to be coupled with 
hf field through the ponderomotive force. The govern- 
ing equations for the nonlinear coupling of these two re- 
sponses will be described later. Furthermore, in Eqs. 
(|TT]) - p^)) . the terms proportional to H are contributions 
from the plasma magnetization current due to intrinsic 
spin of electrons. Notice, however, that when the spin 
effects dominate, the group velocity tends to decrease as 
the frequency of the whistler approaches the cyclotron 
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frequency [20| and/ or k increases its value. The lat- 
ter indicates a negative group dispersion in the propa- 
gation of whistlers. Notice, however, that the nonlinear 
frequency shift precisely depends on the group velocity, 
and can even become larger (when uj approaches ui c ) due 
to the plasma streaming with the flow speed v ez along 
the external magnetic field. 

Let us now compare the two ponderomotive forces, 
which induce slowly varying electrostatic oscillations in 
plasmas. The expressions of these forces, namely CPF 
[21| and SPF [19| [Note that in the expression for the 
spin ponderomotive force given by Eqs. (12) and (13) in 
Ref. [19j , a factor 2 was missing in the denominator, and 
this factor 2 appears in the average of the forces as in 
Eqs. (7) and (8) there for the classical part] acting on an 
individual electron can be written respectively as 



so, rather may be comparable to the spin contribution. 
In this case, ion dynamics might play roles, and also the 
whistler wavelength may be much larger than the plasma 
characteristic length scale. So, this case will not be so 
effective for high energy particle acceleration. We will 
focus on the regimes as discussed above for which the 
excitation of multiple wakehelds is possible. For w > w c 
and x ^ 1) the classical part is also negligible, however, 
this is not relevant to the present study. 

The equations for the If density response which arises 
due to the ponderomotive force of the hf electromagnetic 
field satisfy the electron continuity, momentum balance 
and the Poisson equation, these are respectively 

d t N + 8 z v ez = 0. (16) 
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where £ = z ~ v g t is the co-moving frame of reference 
for the driving whistler. The spin induced ponderomo- 
tive force arises due to the effects of the finite magnetic 
moment of electrons and is obtained by taking the av- 
erage over the fast time scale of the spin force in the 
momentum equation ([2]) [l9| . We note that the spin con- 
tribution to the ponderomotive force is small compared 
to the classical one when the factor \ = hk 2 /m e uj -C 1. In 
this case, one can typically neglect the spin-contribution 
in the linear as well as nonlinear regimes, and thus the 
results will be valid for weakly or strongly magnetized 
low-density plasmas. However, for an exception to this 
rule, see Ref. |19l |. In the opposite limit, i.e., \ H> 1, 
the SPF in Eq. (TT5T) can , indeed, dominate over the CPF 
when the whistler wave frequency is close to the cyclotron 
frequency. In this case, the whistler (with phase speed 
close to c) wavelength (Xw — 2ir/k) can be smaller or 
of the order of the Compton wavelength (Ac = h/mc) 
of electrons, since hk 2 /m e uj scales as \c/\w- This case 
can be relevant in very strongly magnetized (Bo > 10 9 T) 
and very high-density medium (n > 10 36 m~ 3 ) in which 
lu c > ui pe still holds. The whistler wave dispersion can 
then still be linear and the group as well the phase speed 
of the whistlers may approach the speed of light in vac- 
uum. Below, we will see that this is the case in which the 
whistler pulses can maintain their shape over a macro- 
scopic distance, and thus favors the plasma wakefield ac- 
celeration. 

Next, the whistler dispersion may be linear in the 
regimes of Bo ~ 10 s T and no ~ 10 35 m~ 3 where the 
SPF may not be dominant over, but may be comparable 
to the CPF F cz . On the other hand, for ui <C ui c and 
X 3> 1, though the first term in the square brackets in 
Eq. (fT4"|) is negligible, but the second one can not be 
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(18) 



where Ei is the If part of the wave electric field (wake- 
field) and Vf — y/kBTp/rrie is the Fermi speed rel- 
evant for a low-temperature high density plasmas [23| . 

T F = h 2 (37r 2 n ) 2/3 /2k B m e is the Fermi temperature 
with ks denoting the Boltzmann constant. The term 
oc h 2 is the effect of quantum tunneling associated with 
the Bohm de Broglie potential. The ponderomotive force 
contributions are proportional to the constants Ti and T 2 
[H where 
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in which the first terms appear due to CPF [21( and the 
second (oc K) are due to the SPF [19[ • From Eq. (Tl^| we 
note that for x ^ 1, the SPF contribution is substantial 
compared to the CPF, and Fi^ can be approximated as 
oc h that represents purely a spin quantum effect. Equa- 
tions (TT6|) - ([T8l) can be combined to obtain the driven wave 
equation for If perturbations of the Boussinesq-type 
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N = ^d 2 z \E\ 2 - n2d 2 zt \E\ 



where the ponderomotive force contributions are y,\ = 
£o<^p e ri/2TO e o; 2 and [i 2 = eoUJ 2 e kT2/2m e uj 2 in which spin 
effect is hidden. 

Thus, we have a new set of three coupled equations, 
namely dTUJ) , CU) and $ZU§, modified by the SPF and 
the effect of quantum tunneling, which describes the 
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nonlinear coupling of electron whistlers with the field 
aligned electrostatic density fluctuations. These equa- 
tions can be written in the nondimensional forms as 
(z -> z/X Fe ,t -> tu> pe ,E -> J5/£b, u ez -> v ez /c s , where 
Afc = Vp/uJpe = X e V F /c) 
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and, in particular, in absence of the quantum force, 
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[d 2 t + H 2 di -d 2 z + l]N = \E\ 2 - C 2 dl\Ef, (23) where K = 1/{V 2 - 1)W 



where £7 = ^2k B T F n /e , V g 



v g /V F ,V' 



v'gUJpe/Vf, V T / = A/o; pe , H = hujp e /2kBTF is the quan- 
tum parameter and £i = ^peTi/u^ 2 , £2 = WpgfcVF^/w 2 



Equations (f2Tj) - (|23|) contain the main results of the 
present work. 



III. EXCITATION OF THE WAKEFIELDS 

In this section, we consider the excitation of one- 
dimensional wakefield by the hf whistlers propagating 
with the group velocity. Thus, we look for station- 
ary solutions of our main Eqs. (|2~Tj) - (|23|) in the frame 
£ = z — V g t. To this end, we assume E to be of the form 
E = W (£) exp (—iOt) , where W is a real function nor- 
malized by Eq and is a real constant. Then Eqs. (f2Tj) - 
(l23l) reduce to 



d 2 W 



WQ - ANW = 0, 
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[H 2 d\ + (V g 2 - 1) rf 2 + 1] N = (Ci + ( 2 V g ) d\W\ (25) 



where 
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(26) 

Equation (1231) shows that the electrostatic wake field is 
created by the self-consistent whistler waves. We also 
find that the inclusion of the quantum correction cx H 
leads to two different characteristic length scales, which, 
in turn, will give rise multiple oscillatory wakefields. 
Thus, considering a trial solution of the corresponding 
homogeneous equation of (|25[) as N ~ exp(ifc p £), we ob- 
tain kl - n 2 a k 2 + 4 = 0, where K a = (V 2 - if^/H 

and «b = 1/viF, This biquadratic equation has, in gen- 
eral, two pairs of real roots ifcj.,2 *2] for = v g /Vp > 
y/i + 2H where 



1.2 — _2. x a / 
K l,2 - o + V , 



(27) 



For a Gaussian driving pulse of the form W ~ 
exp (— £ 2 /L 2 ) , Eq. (|2"5j) yields the wakefield behind the 
driving pulse, i.e., for |£| ^> L p as AT = N% — Ni where 



Next, for the generation of whistler driven wakefields, 
we numerically solve the Eqs. ((24)) and (1251) by Newton's 
method with the boundary conditions AT, W, dN/d£, 
dW/dt; — > as £ — > ±00. In order to study the influence 
of the quantum tunneling effect, we consider two cases, 
namely for H = [Note, however, that H (~ 1/K) = 
does not mean that one recovers the h = case, rather 
it implies that we have simply disregarded the quantum 
tunneling effect associated with the Bohm potential] and 
for nonzero H. Here the typical length scale for plasma 
collective oscillation is the Fermi screening length Xp e = 
X e V F /c < X e as long as the Vp < c, and it must be such 
that Xfc < Xw, the whistler wavelength. The typical 
time scale is the electron plasma period, which is shorter 
than the ion plasma period for which the ion dynamics 
is negligible. Moreover, quantum effects become impor- 
tant when C = T F /T c i asslca i = (1/2) (3^ 2 n A|) 2/3 > 1, 
where Xb = h/m e V F is the thermal de Broglie wave- 
length. Furthermore, the quantum coupling parameter 

9q ~ (l/ ri oA'p e ) 2 ^ 3 < 1 defines a quantum collisionless 
regime where the collective and mean-field effects domi- 
nate. We now study the two cases as follows. 

In the case of H = 0, the plasma characteristic length 
scale for the wake is X p = 2it/k p , and we choose the pa- 
rameters as B a = 3 x 10 8 T, n a = 10 35 m~ 3 and u = 0.78, 
= 0.22, for which lj/uj p — 2.3, u c /ujp — 3,v g /Vp = 
l.3,v g /v p = 0.7, X Fe /X e = 0.4 and E = 1.7 x 10 16 V/m. 
In this case the plasma oscillation length is A p (= 5.4) > 
X\y (= 2n/kX Fe in normalized form) = 5.3, the group 
dispersion V g = —2.7 and the nonlinear frequency shift 

A = 0.4. Moreover, g Q = 0.11 < 1 and C ~ 1- Thus, 
in the quantum regime, the short-scale wakefield can be 
generated as shown in the upper panel of Fig. 1. In- 
creasing the magnetic field strength, namely Bq — 6 x 10 8 
T and keeping electron concentration as the same, i.e., 
no = 10 35 m~ 3 , we find that the group velocity increases 
(v g /vp = 0.85) and hence for a increasing value of the 
ratio u c /ujp — 5.9, the wakefield remains coherent, but 
its amplitude becomes much lowered (not shown in the 
figure). The whistler group dispersion has now been 
lowered as V' = —1.5 than the previous case and the fre- 
quency is shiftcd-up to A = 0.53. The other parameter 
values in this regime are, e.g., X p — 10.9 and Xw = 3.2. 

The lower panel of Fig. 1 shows the wakefields for a 
nonzero H. We find that the plasma wake electric field 
is enhanced due to the quantum tunneling effect. Here 
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FIG. 1: (Color online) Plasma wakefields [Numerical solution 
of Eqs. (|24[) and (]25p] driven by the whistlers (thick line) for 
H = (upper panel) and for H — 0.13 (lower panel) with 
characteristic length scales X p = 5.4 (A.Fe) and Ai = 5.5 
(Af e ) respectively. The case in which SPF is comparable to 
CPF. The parameter values are given in the text. 



the corresponding plasma characteristic length scales are 
Ai,2 = 2tt/ ki t 2 as defined above. For the length scale 
Ai, we use the parameters B = 6 x 10 8 T, no = 2 x 10 35 
in" 3 , uj = 0.78 and 6 = 0.04 for which uj/uj p = 3.3, 
uj c /uj p = 4.2, v g /V F = 1.3, Vg/vp = 0.77, k\ Fe = 1.9, 
H = 0.13, Ai = 5.5, X w = 3.3, V' g = -1.5 and A = 0.5. 
We note that the length scale A2 docs not favor the 
wakefield generation in the parameter regimes considered 
above, since for V g > y/T+~2H, A2 scales as the Compton 
wavelength, which becomes lower than the whistler wave- 
length Xw- However, since as said before, Xw may ap- 
proach or lower than the Compton wavelength in the case 
of very strongly magnetized and very high density plas- 
mas where the SPF dominates over the CPF, multiple 
wakefield generation can be possible as can be illustrated 
from Fig. 2 below. Thus, for the parameters Bo = 4x 10 9 
T, n = 10 36 itT 3 and ui = 0.78, 8 = 0.008 for which 
uj/uip = 9.7, uj c /(jjp — 12.5,v 9 /Vf = l-l,w 9 /w p = 0.97, 
kX e = 9.97, kX Fe = 8.4, 27r/fcA e = 0.63, H = 0.13, 
Xw — 0.75, Vi — —0.1, A = 0.3, the corresponding wake- 
field excitation is shown in the upper panel of Fig. 2 with 
Ai = 2.8. From the lower panel of this figure one can see 
the similar excitation but with a shorter scale A2 = 1.45. 
In the latter, wakefield amplitude is seen to be enhanced 
compared to the case with Ai. 
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FIG. 2: (Color online) Plasma wakefields [Numerical solution 
of Eqs. (124[) and (] 25 P ] driven by the whistlers (thick line) for 
a fixed H = 0.13 but for two different characteristic length 
scales Ai = 2.8 (Af c ) (upper panel) and A2 = 1.45 (A_p e ) 
(lower panel). The case in which SPF dominates over CPF. 
The parameter values are given in the text. 



IV. DISCUSSION AND CONCLUSION 

Numerical results in the previous section suggest that 
wakefield excitation by the whistlers can be possible in 
a strongly magnetized dense plasmas ( Bq < 10 8 T and 
n o 10 35 m~ 3 ) with intrinsic spin of electrons. In this 
regime the SPF is no longer negligible rather may be 
comparable to the CPF. Such parameter regimes can be 
achievable in the interior of white dwarf stars where elec- 
trons are indeed degenerate and the electron degeneracy 
pressure is what supports a white dwarf against gravita- 
tional collapse. The magnetic fields in such compact de- 
generate stars might be due to conservation of total sur- 
face magnetic flux during the evolution or by the emis- 
sion of circularly polarized light or else. However, a small 
number of white dwarfs have been examined for magnetic 
field strength exceeding 10 5 T. 

In absence of the quantum force, the wakefields can 
be generated with the characteristic length scale X p = 
2w(Vg ~ l) 1 / 2 and in the parameter regimes Bo ~ 10 8 
T, n ~ 10 35 m -3 as described above. Apart from the 
Fermi pressure, the quantum tunneling effect associated 
with the Bohm potential introduces an additional higher- 
order dispersion. The latter, in turn, introduces two char- 
acteristic length scales Ai and A 2 (< Ai) (quite distinctive 
from the case with H = 0) out of which only Ai favors 
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the wakefield excitation in the said regimes, since A2 be- 
comes smaller than the whistler wavelength \w (> 10~ 12 
m), a case not effective for the wakefield generation. 
However, multiple wakefield excitation can be possible 
with two length scales Ai and A2 in very strongly mag- 
netized and superdense plasmas with Bq > 10 9 T and 
n o ^ 10 36 ni -3 . In these regimes the CPF is strongly 
dominated by the SPF. The latter will then dominantly 
accelerate the electrons by separating the electric charges 
and building up a high electric field. Furthermore, such 
parameter regimes could be relevant in the deeper layer 
(above the neutron-drip layer) of the crust of neutron star 
in which electrons are degenerate implying that electri- 
cal and thermal conductivity may be huge because the 
electrons can travel great distances before interacting. 
However, the physics and the composition of neutron 
star interiors are not yet fully understood. 

We ought to mention that in the regimes of very strong 
magnetic fields and superdense plasmas as mentioned 
above, the nonrelativistic fluid model may no longer be 
appropriate as in such cases the Fermi speed approaches 
(or can be greater than) the speed of light in vacuum and 
whistler wavelength may become smaller than the Comp- 
ton wavelength of electrons. In this situation, relativis- 
tic spin-quantum fluid model or kinetic approach could 
be well set. By increasing the magnetic field strength 
or the ratio lu c /lo p , we see that whistler wave can have 
smaller group dispersion. This could be important for 
whistler waves to be a viable mechanism in magnetized 
dense plasmas in order to gain higher energy for an ac- 
celerated particle from the plasma wakefield. 

Next, as an alternative mechanism instead of using 
lasers or particle beams as drivers which have been used 
mostly for laboratory plasma wakefield accelerations, the 
EM whistler waves should be of fundamental interest in 
plasmas, in particular, quantum plasmas for astrophys- 
ical settings. However, conclusive evidence needs fur- 
ther investigation in this area by considering, e.g., full 
scale simulation of a relativistic fluid model or kinetic 
one accounting for the quantum statistical and mechani- 
cal (tunneling) effects as well as intrinsic spin of electrons. 
The latter might have a role in the regimes previously 
considered to be as classical [3l[. We mention that our 
model can be well applicable for low or moderate density 



magnetized plasmas as well where electrons may not be 
degenerate, spin effects might not be so important. In 
this case one can consider, e.g., the isothermal equation 
of state instead of the Fermi pressure law for electrons. 
Thus, classical results can also be recovered by disregard- 
ing the term oc h for the spin contribution and the term 
oc H 2 for the quantum tunneling effect. 

To mention, neutron stars are known to be compact 
and carry intense surface magnetic fields about 10 9 T or 
more .32]. It has been investigated that short gamma- 
ray bursts (GRBs) may arise (though, there are still a 
lot of uncertainties about their origin) from collisions be- 
tween a black hole and a neutron star or between two 
neutron stars. However, when neutron stars collide, the 
tremendous release of energy results into highly relativis- 
tic out-bursting fireballs (jets) [Hj]. The latter are most 
likely in the form of a plasma and can have initial plasma 
density > 10 32 m -3 . Such collisions of intense magnetic 
fields may create strong magnetoshocks where whistler 
waves are embedded. 

To conclude, we have presented a theoretical investi- 
gation for the possible excitation of wakefields by self- 
consistent whistler wave field at nanoscale in a magne- 
tized spin quantum plasma. The whistler wave enve- 
lope is governed by a NLS equation coupled to a driven 
Boussinesq-like equation for the If wakefield. The quan- 
tum force is shown to be responsible for the excitation 
of multiple wakefields in strongly magnetized superdense 
plasmas in which SPF strongly dominates over the CPF. 
Furthermore, the effect of quantum tunneling is to en- 
hance the wakefield amplitude, however it is reduced by 
the external magnetic field. Finally, the present investi- 
gation can be generalized to multi-dimensional wakefield 
excitation in relativistic spin quantum plasmas. 
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